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Abstract 



A possibility to search for T-violation in neutrino oscillation with the use of 
k> , muon polarization is studied. The sensitivity to T-violation is examined with various 

H I magnitudes of muon polarization as a function of muon energy and long-baseline 

distances. 

1 Introduction 

In this note, we discuss a possibility to discover the CP-violating effect in the leptonic 
sector [H] at a neutrino factory. Observation of CP- violation would imply the measurement 
of an imaginary part of the couplings in the Lagrangian associated with leptons 0. There 
are two ways to study CP violation. They are 

1. to study a difference between the CP conjugate modes, or 

2. to study a difference between the T conjugate modes. 
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The first method is to observe the difference between the particle and its anti-particle. 
In a neutrino factory, it can be studied by comparing the difference of the appearance 
event rates between z/e -^ v^ and z/g -^ t'^ [S §; i, li 0; H 3 • However, it has been pointed 



out that the matter effect would introduce a sizeable fake CP-odd effect 10, O . It is 



necessary to discriminate the genuine CP violation effect from the fake matter effect. 

The second method is to observe the difference between the transitions v^ — > z/^ and 
i^fj. ^ ^e- By the CPT theorem, in vacuum, observation of T- violation is identical to 
observation of CP-violation. Even in a terrestrial oscillation experiment, since a possible 
T-odd effect from asymmetric matter density profile in the earth is very small [1^, this 
T-reversed difference reveals as a very clean signal B O] on the imaginary phase in the 



lepton sector R, |, 10, 12, 13, 14, 15, 0, 17, 18, 19, 20, 21, 22, 53, 24, 25, E 



Therefore, the comparison of T-reversed oscillation modes would give clear signature 
than those of CP-reversed oscillation modes. However, although it is easy to observe 
t'e —>■ ^n oscillation, it is known to be difficult to observe v^ — > z/g oscillation at a neutrino 
factory. The reason is as follows. At a neutrino factory, v^ is generated by the decay 
of muons as well as Ve- The former after the oscillation from v^ -^ z/g, produces e~s 
and the latter from muon decays produces e"^s at a detector. It is difficult to distinguish 
e^ from e~ after the creation of electromagnetic shower in most of the detectors being 
considered at present. Therefore, it is difficult to identify the oscillation events from the 
non-oscillation events which is larger in statistics. 

It is critical to examine whether there are any solutions to solve this problem for 
the search for T-violation modes useful. There was one idea to identify the oscillation 
events by the use of polarized muons in the muon storage ring ||T^. The principle of the 
idea is the following. Supposed that the muons in the muon storage ring have their spin 
polarization P^- = — 1, there would be no Ue and only z/^ at the very forward direction 
along the muon momentum direction. Observation of e— like events at a detector would 
directly imply the appearance of u^ -^ z/g. In reality, the muon spin polarization would 
not be 100 %. But, if the spin polarization of the muons in the muon storage ring can be 
changed, accordingly the neutrino energy spectrum shapes can be changed. The changes 
of the spectrum of neutrino and anti- neutrino are different. Therefore, in principle we 
could discriminate the events associated with neutrinos from those from anti-neutrino. In 
turn, the oscillating events and non-oscillating events can be separated. 

2 Principle 

In this section, the principle to discriminate the oscillating events by using the muon 
polarization is presented. The neutrino flux in the very forward direction of a fi~ beam 
with its muon spin polarization P^ is given by 

$p, oc y\l-y)-{l + P,), (1) 

$,„ oc y\3-2y) + P^-y\l-2y), (2) 



where y = Ei,/E^. E^ is the fi~ energy and E^, is the neutrino energy. For /i"*", eqs.([5) 
and (0) can be given with changing the sign of P^. From using eqs.(|l|) and (0), the 
non-oscillating event rate of z/g — > i^e {Np^^pJ and the appearance event rate of u^ -^ v^ 
{N^.-.uJ are given by 



iVp^^p, = (l + P^)iV3, (3) 

N,^^,^ = N, + P,N,, (4) 

respectively. Here A'^i, N2 and A^3 are the event numbers corresponding to each term of 
neutrino flux. They are typically represented by 

E^ /-SmaxZ-B^ 

Ni = C-^ dy2y\3-2y)P,^^,SE,)a{E,), (5) 

N2 = C^ r^^^^^""" dy2y\l-2y)P,^^,XEu)(r{E,). (6) 

^3 = C^ r"^^^""" dyl2y\l-y)P,^^pXE.)cT{E,)^ (7) 

where P^^^^^, Pp^^Pe are oscillation probabilities, a (a) is the detection rate including 
the neutrino (anti-neutrino) cross section and signal efficiency, Ej^m (-E'max) is minimum 
(maximum) neutrino energy observed in this assumption, L is the baseline length, and C 
is the constant proportionate to the decay muon number and the detector mass. From 
these equations, the rate of e— like events, which is the sum of Np^^p^ and N^^^i^^, can be 
given by 

iVe-like = /3o + ^^/3l, (8) 

/3o = Ni + Ns, (9) 

A = N2 + N,. (10) 

Note that for P^ = —1, A^e-uke becomes equal to Ny^^y^ itself from eq.(^. The estimate 
of the e— like event rate at P^ = —1 is essentially the measurement of N^^^y^. In another 
words, if (3o and /5i can be estimated, the direct observation of Ny^^y^ could be made in 
principle. This is shown in Fig.|l]. 

3 Procedure 

In this section, the procedure to estimate the appearance event rate, Ny^^y^, at P^ = —1 
from the measurements is given. In reality, the muons in the muon storage ring have 
some distribution of their polarization. It is assumed that a number of the muons of their 
muon polarization of P^^, {i = 1 ■■■ n) is Ni{= N^/ f-^, {i = 1- ■ -n). They yield yi e— like 



K 



e— like 




Figure 1: Graphical view of the idea. We have observations of e— like event with partially 
polarized muon beams and from them we estimate the event rate at P^ = — 1. 

events. A^^ is a total number of the muons in the ring. For those, eq.@ can be rewritten 
by 



(11) 
(12) 

-1-P,„ (13) 

(14) 
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(15) 



The expectation values of yt, Hi is given byQ 



f^yi = f3o + P,*Pi- (16) 

It is assumed that ?/j's are large enough so that it would follow the normal distribution 



^In the following, we denote /?q by Pq. 



of N(?/j, Hi). In this case, the hkehhood function to estimate Po and Pi is given by00 

^ {/.w-(ft + P,.-ft)}' 

-'°«^=? — ^. — ■ <"* 

From this hkehhood function, the equations for the estimates of Pq and Pi can be give by 



^.: , o ^ (^/ 



*\2 



/^oEt^t + z^iE 



A«j. 



Jj 2/i j Ji Hi 

Here the definitions of P^* by eg. (]TB|) , and of P^ by eq.(|T^ are used. Also the last 
term in the first line of eq.(|19D is removed. Using these eqs. (p!8D and (0), Po and Pi in 
terms of i/i can be estimated as follows: 

F (P *)'^ 
/5o = nE^f^ and (20) 



Dr ftvi 



F P * 

/^i = -7^E7£7' (21) 



where 



~\rfh^)[r fh^ ) [rfhj 



(22) 



Thus, Ny^^y^ at P^ = — 1 can be estimated using eqs.(pU[) and (^) as 

K^_,,^ = P, + {-l-P,)Pi. (23) 

The statistical error on this estimated value is now evaluated. Formally, the variance 
of the estimate is given by 



^Exactly speaking, the denominator of eq.(p7|) ffvi- However, for simplicity, under the assumption 
that t/i's are large, it can be replaced with ffyi. This is justified since the non-oscillating event rates 
Njy^^jj^ at any polarization are large. (It does not depend on whether there are enough appearance event 
rate N^, ^^^ as long as the polarization is not very close to —1.) 

■^The likelihood function with Poisson distribution can be also treated. By this, the same function as 
eq.(O) can be obtained under the assumption that j/^'s are large enough. 



V{N. 



V{(3o + (-1 - P,)(3i) = V{(3,) + (-1 - P^) V(A) + 2(-l - P^)V{f5o, A),(24) 



where V{[3) is the variance of /3, and V^(/3o,/3i) is the covariance of /3o, A- 

It is impossible to evaluate the variance of N^^^^e using eqs. pO|) and ( pT]) directly. 
Therefore, it is approximated by the inverse matrix of the expectation value of Fisher's 
information.^ Fisher's information matrix is given by 



92 



'e^ 



dPidp, 



logL 



fhi 



and hence the variance matrix is given by 



y^ ^H 

i Ji Vi 
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I (P *V 

}, Ji Hi 
p * 



P * \ 



i Ji Hi 
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ffVi 
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fiVi J 



Thus, the estimate of N^^^y^ has its variance of 

1 



\ 



flVt ) 



v{N,,^,:) 



D\^ fh^ 



+ 




1 



2N„ 



+ Nf 



-1 - Pn 



.2 E.i//^. 



Y.^{P,:Y I ffy^ 

Here Np is the estimate of the event rate at P^ = P^ when A*"^- /i~s are used. 



(25) 



(26) 



(27) 
(28) 



4 Numerical Analysis 

In the analytical treatment shown above, several assumptions are made. For instance, it 
is assumed that the likelihood function of the /5's is given by eg. (pTj) , and the variance of 
the estimate is given by the inverse of the Fisher's information matrix. 

In this section, numerical calculations is given. What we have to calculate numeri- 
cally is the extrapolated value at P^ = — 1 which indeed coincides with the theoretical 
expectation value y\p^^_i and the variance of it which is actually given by eg. (pS]) . 

The algorithm for this numerical analysis is as follows: 



''strictly speaking, the inverse matrix of the expectation value of Fisher's information gives the lower 
limit of the corresponding variance. However if the likelihood function is constructed well then it coincide 
with the variance and hence we use this approximation. We also approximate the expectation value of 
the Fisher's information by eg. (Pq) since it is assumed that yiS are large number enough. 
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step 0. Select several values of the polarization P^^ for the measurements, 

step 1. Fix all the parameters, such as the theoretical parameters (e.g. mixing angle), a 
muon energy and so on, and then calculate the expectation of event rates in each 
polarization values (|/j). 

step 2. Generate the "actual" events of yi according to Poisson distribution with its expec- 
tation value yi 



step 3. Substitute these event sets into eqs.pO|), pT|), and find the extrapolated value at 
P^ = —1 from eq.(|23|) by using these /5's. 



step 4. Iterate the step 2 and step 3 in several times to find the type of distribution of the 
extrapolated value. 

This virtual measurements for various parameters were made. It is found that the extrap- 
olated value follows the normal distribution with its mean value y\p =_i and its variance 
given by eq.(|28|). 



5 Sensitivity for T Violation 



In this section, the sensitivity to T-violation using the present procedure is discussed to 
find the statistics to see the T-violation effect. 






where 



N : AT,^^,^, 

and the subscript indicates the estimate of the event rate with CP violating phase 
S = 5o = {0,7f} and subscript i represents ith energy bin. Note that in the denominator 
of x^ there is V. Although V is given by eq.(p8D theoretically, it is calculated by using 
the numerical method in section ^. To confirm that T violation can be observed at 90% 
confidence level, x^{^) > X9o%('^) i^ required. Suppose that three measurements can be 
adopted with different muon polarization values (P^^ = {— P^,0, P^}), the sensitivity of 
90% confidence level can be calculated. It is shown in Fig.^ We chose the most sensitive 
method among n = 1, 3, 5 in each {E^, L) region. The usage of statistics is referred in 
Ref.||ll|| in detail. 



Here, it is discussed how many neutrinos from muon decays should be needed to 
observe the T-violation effect. First of all, for simplicity, iVj = N^+ is assumed. Here, 
the following theoretical parameters are used; sin 6^12 = 0.5, sin 6^23 = 1/v^, sin 6^13 = 0.1, 
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Figure 2: Required data sizes in the unit of 10^^ niuon decays x 100-kt detector to obtain 
a sensitivity of 90% confidence level are plotted as a function of muon energy and baseline 
length for various sets of muon polarizations P^^. The meaning of each line is explained 
in Fig.^. For example, the contour with the label of 0.5 implies 5 x 10^° muon decays with 
a 100-kt detector is used. For each case, three measurements with different Pa values 
are assumed. These three values of the muon polarization are shown for each plot. If the 
muon polarization is high, the necessary data size is smaller. 
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Figure 3: Same as Fig.0, but here bm^^ = 5 x 10 eV 



5o\^2 



5mli = 3 X 10 ^eV^, 5mli = 10 ''eV^, 5 = 7r/2. We assume that each parameter including 
matter density has 10% ambiguity. The sensitivity is almost proportional to 

sin 5 sin 2^i2 sin 2^23 cos ^13 -— |i . (30) 

Therefore, the sensitivity becomes 4 times worse when Sm^i = 5 x 10~^eV^ (see FigJ^). 

In Fig.0, it is found that at worst P^ = ±0.3 should be needed to observe T-violation 
effect with a total of 10^^ muon decays and a 100-kt detector. On the contrary to the sen- 
sitivity of CP-violation effect usually discussed, the sensitivity will not change drastically 
by including other sources of errors, since the V will convey the largest error among the 
other uncertainties. 

In the remaining of this section, we see what we can learn from the theoretical analysis. 

5.1 How many fi'^ and ^u decays are needed ? 

In this subsection, the most efficient ratio, fi, between the numbers of fi~^ (^^+) and /i~ 
(A*"^- ) to obtain the best sensitivity to CP violation is estimated for a fixed distribution 
of A*"^- at each polarization P^^. To do this, the x^ can be rewritten to factor out the 
numbers of events as follows, 



N ^ N^+N, 

(NqN - NpNf 
^ ^ N^V/N,-+N^N/N^^- ^ ^ 

If the total numbers of muon is given, then 

A^ - 
J2Ni + N^+=Y,^ + N^+= FN^- +N^+=N (32) 

is the given number of the total events. Here, F is given by eq.(|l^). Under these condi- 
tions, the most efficient ratio of the numbers of fi'^{N^+) and fi~{N^-) is given by 



N^^ ---./- (33) 



where A is given by eq.(^) xF 
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Figure 4: Uniform Distribution of polarization of muon, being distributed from — P^ to 
+Pfj, uniformly. 

Since the oscillation probability is at most a few %, while Np^ in eq.(^) is essen- 
tially given by the non-oscillation event, from the fact that N/Np ~ O(10~^) and that 
J2i[^/{ffyi)]/J2i[{PiJ.*i)'^/{fiyi)] in eg. (P^) is much larger than a unity in a realistic case, 
it is found that a much larger number of /i~ is needed. Therefore by tuning the ratio of 
the numbers of /x"*" and /i~, a higher sensitivity to T- violation effect can be obtained by 
several tens %. 

5.2 Dependence of Muon Polarization Distribution 

The measurements of Ni events with muon polarization P^^ are usually made. In this 
subsection, the optimization of A^jS {i = 1 ■ ■ -n) for each given P^^ to obtain the best sen- 
sitivity is studied. This is equivalent to examine the best muon polarization distribution. 

For the most simple case to consider is a uniform distribution of Ni from — P^ to +P^ 
as shown in Fig.^. 

Let's consider two cases. One is a uniform distribution of N^ from — P^ to +P^ as 
shown in Fig.^ The other is discrete measurement at ±P^ as seen in Fig.^ 

It is assumed that the same numbers of muon decays are used for the two different 
distributions of Fig.§ and Fig.^ In this case, the variances (p8|) of the extrapolated 
appearance numbers are different, depending on the distribution. In case of Fig.^,Q 

^Here we assume that we do not have very large polarization. It means that j/i's are almost constant. 
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Figure 5: Distribution of the polarization of /j,. There are only muons with their polar- 



ization around — P^ and +P^. 
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^^^/^^^^ loc-^-1, (34) 



i^E.(p.:)7/?i/. p, 



fj- 



while in case of Fig.^ 



Eil//?l/i _io, 1 _i_ (35) 



Therefore, a three times higher sensitivity can be obtained for the distribution of Fig.^. 
In reahty, there will be a loss of muon decays to create the distribution in Fig.^. If this 
reduction factor of the number of events is defined to be M (namely only (1/M) events 
can be used), the sensitivity is given by 

Then the sensitivity is better by 3/M. To discuss further, it is needed to have realistic 
muon polarization distribution. However, once it is obtained, the discussion to study a 
sensitivity is quite straightforward. 

6 Summary and Discussion 

In this note, we have examined a possibility to study T-violation with the use of muon 
polarization. The physics motivation to search for T-violation is to observe direct CP 
violation arising from the imaginary phase in the neutrino mixing matrix. The advantage 
of search for T-violation is to have a negligible small matter effect. 

In the present study, the sensitivity to T-violation has been studied as a function 
of muon energy and long baseline length for various set of muon polarizations. In each 
set, it is assumed that at minimum three measurements with different values of muon 
polarization are used, for instance, P^^ = ±P^ and are used in the present study, where 
P^ is the magnitude of the muon polarization. Needless to say, it is shown that the muon 
polarization achievable is larger, the sensitivity is better. From the present study, it is 
found that |P^| > 0.3 is necessary when it is assumed that we can get a total of 10^^ muon 
decays and a 100-kt detector. 

The magnitude of muon polarization is being studied in the accelerator study. It is not 
easy to have large muon polarization in a very high intensity. It would be needed that the 
magnitude of muon polarization and a number of muons have to be traded off. However, 
if search for T violation can be made with realistic magnitude of muon polarization, it 
would provide unique way to study the CP phase in the MNS leptonic mixing matrix. 
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